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SUMMARY 


High performance elements are simple finite elements constructed to deliver engineering accu- 
racy with coarse arbitrary grids. This paper is part of a series on the variational basis of high- 
performance elements, with emphasis on those constructed with the free formulation (FI 1 ) and 
assumed natural strain (ANS) methods. The present paper studies parametrized variational prin- 
ciples that provide a foundation for the FF and ANS methods, as well as for a combination of 
both. 


1. INTRODUCTION 

For 25 years researchers have tried to construct “best” finite element models for problems 
in structural mechanics. The quest appeared to be nearly over in the late 1960s when 
higher order displacement elements dominated the headlines. But these elements did 
not dominate the marketplace. The overwhelming preference of finite element code users 
has been for simple elements that deliver engineering accuracy with coarse meshes. The 
search for these “high-performance” (HP) elements began in the early 1970s and by now it 
represents an important area of finite element research in solid and structural mechanics. 
Many ingenious schemes have been tried: reduced and selective integration, incompatible 
modes, mixed and hybrid formulations, stress and strain projections, the free formulation 
(FF), and the assumed natural strain (ANS) method. 

The present paper is part of a series [8-12) that studies how several high performance 
element construction methods can be embedded within an extended variational framework 
that uses parametrized hybrid functionals. The general plan of attack is sketched in Figure 
1. Heavy line boxes are those emphasized in the present paper. The extensions, shown on 
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the left, involve parametrization of the conventional elasticity functionals and treatment 
of element interfaces through generalizations of the hybrid approach of Pian [14-16]. 

The effective construction of HP elements relies on devices, sometimes derisively called 
“tricks" or “variational crimes,” that do not fit a priori in the classical variational frame- 
work. The range of tricks range from innocuous collocation and finite difference constraints 
to more drastic remedies such as selective integration. Despite their unconventional na- 
ture, tricks are an essential part of the construction of high-performance elements. They 
collectively represent a fun-and-games ingredient that keeps the derivation of HP finite 
elements as a surprisingly enjoyable task. 

The present treatment “decriminalizes” kinematic constraint tricks by adjoining La- 
grange multipliers, hence placing the ensemble in a proper variational setting. Placing 
formulations within a variational framework has the great advantage of supplying the gen- 
eral structure of the matrices and forcing vectors of high performance elements, and of 
allowing a systematic derivation of classes of elements by an array of powerful techniques. 

Note the reliance of the program of Figure 1 on hybrid functionals. The original 
1964 vision of Pian [14] is thus seen to acquire a momentous significance. It is perhaps 
appropriate to quote here the prediction of another great contributor to finite elements: 


T. H. H. Plan responded to the problem of plate bending by Inventing the 
"hybrid formulation", which avoids the problem of slope continuity. He 
assumed that the element responds not according to shape functions but 
according to element stress fields. These communicate with the outside 
world via the boundaries .... Hybrid elements can be the most competitive 
and we believe that the future lie In that direction. However, the formula- 
tion Is more complicated. Therefore we advocate that researchers should 
try to cajole their formulation Into shape function form, so that users do 
not have to struggle. In the form, hybrid elements are no more difficult 
to use than the Iso-P elements ... Unfortunately at the time of writing 
we have no uniform technique to achieve this. 

B. Irons and S. Ahmad, Techniques of Finite Elements (1980), p. 159 


Fulfillment of the prophecy appears to be near. 

2. THE ELASTICITY PROBLEM 

Consider a linearly elastic body under static loading that occupies the volume V . The 
body is bounded by the surface 5, which is decomposed into S : Sd U St. Displacements 
are prescribed on Sj whereas surface tractions are prescribed on S t . The outward unit 
normal on S is denoted by n = n,-. 

The three unknown volume fields are displacements u = infinitesimal strains e = e*,-, 
and stresses a = a»y. The problem data include: the body force field b = 6, in V, prescribed 
displacements d on S.j, and prescribed surface tractions t = on St. 
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Projections 


Collocation 


Lattice Treatuent 



Figure 1 Program of attack for variational formulation of HP elements 


The relations between the volume fields are the strain-displacement equations 
e = ^(Vu + V T u) = Du or e,-y = + «y ,i) in V, 

the constitutive equations 

o E e or 0»y 
and the equilibrium (balance) equations 

— div a = D*<t = b or aijj + = 0 in V, 
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in which D* = -div denotes the adjoint operator of D = j(V + V r ). 

The stress vector with respect to a direction defined by the unit vector v is denoted as 
0 ^ — €r.v, or Oyt — (TjjVy. On S the surface-traction stress vector is defined as 

a n = <r.n, or = Oijtij. (4) 

With this definition the traction boundary conditions may be stated as 

a n = t or °ij n i — U on ( 5 ) 

and the displacement boundary conditions as 

u = d or Ui = di on Sj. (6) 

3. NOTATION 


8.1 Field Dependency 

In variational methods of approximation we do not work of course with the exact fields 
that satisfy the governing equations (1-3, 5-6), but with independent (primary) fields, which 
are subject to variations, and dependent (secondary, associated, derived) fields, which are 
not. The approximation is determined by taking variations with respect to the independent 

fields. 

An independently varied field will be identified by a superposed tilde, for example u. 
A dependent field is identified by writing the independent field symbol as superscript. For 
example, if the displacements are independently varied, the derived strain and stress fields 

are 

e « = i(V + V T )u = Du, ff “ = Ee“ = EDu. (7) 

An advantage of this convention is that u, e and a may be reserved for the exact fields. 

8.2 Integral Abbreviations 

Volume and surface integrals will be abbreviated by placing domain-subscripted paren- 
theses and square brackets, respectively, around the integrand. For example: 

V)yVj v fiV, \f)s=J s fdS, fjdS, I f\s'=J s / dS - <»> 

If f and g are vector functions, and p and q tensor functions, their inner product over V 
is denoted in the usual manner 

(f,gV = f / f g dV= [ fi 9 id.V , ( P ,q)v = f / P <1^ = / PijQij^ (») 

Jv Jv Jv 

and similarly for surface integrals, in which case square brackets are used. 
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Figure 2. Internal interface example. 

8.8 Domain Assertions 
The notation 

(a = 6)v, [a = 6]s, [a = *1^, [<* = &]s«, ( 10 ) 

is used to assert that the relation a = b is valid at each point of V, S, S,i and S ti respectively. 
8.4 Internal Interfaces 

In the following subsections we construct hybrid variational principles in which boundary 
displacements d can be varied independently from the internal displacements u. These 
displacements play the role of Lagrange multipliers that relax internal displacement con- 
tinuity. Variational principles containing P d will be called displacement-generalized , or 
d-generalized for short. 

The choic e of d as independent field is not variationally admissible on S 4 or S t . We 
must therefore extend the definition of boundary to include internal interfaces collectively 
designated as Si. Thus 

S : S ti U S t U S { . (11) 

On Si neither displacements nor tractions are prescribed. A simple case is illustrated 
in Figure 2, in which the interface S’,- divides V into two subvolumes: V + and V~ . An 
interface such as Si on Figure 2 has two “sides” called S+ and Sf , which identify Si 
viewed as boundary of V + and V ~ , respectively. At smooth points of Si the unit normals 
n + and n“ point in opposite directions. 

The integral abbreviations (8)-(9) generalize as follows, using Figure 2 for definiteness. 
A volume integral is the sum of integrals over the subvolumes: 

(f)v =' \ fdV+f fdV. (12) 

Jv+ Jv- 
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An integral over 5 ,- includes two contributions: 



where g + and g~ denotes the value of the integrand g on S+ and S i , respectively. These 
two values may be different if g is discontinuous or involves a projection on the normals. 

Following a finite element discretization, the union of interelement boundaries becomes 

4. THE ELASTICITY FUNCTIONALS 
The variational principles of linear elasticity are based on functionals of the form 

n = U - P, (14) 

where U characterizes the internal energy stored in the body volume and P includes other 
contributions such as work of applied loads and energy stored on internal interfaces. We 
shall call U the generalized strain energy and P the forcing potential. 

It must be pointed out that all functionals considered here include independently varied 
displacements. Thus, the class of dual functionals such as the complementary energy are 
not included in the following study. 

4.1 Volume Integrals 

The generalized strain energy has the following structure: 


U = \jix{o,n a )v + ji 2 &,e)v + ii3(<7> e “)v + |j22(^ e ,e)v + j 23 (*‘\e“)i' + 5.733 (* M ,e u )v 

(15) 


where ju through j '33 are numerical coefficients. For example, the Hu-Washizu principle 
is obtained by setting J12 = — 1 > Jia = 1 » J22 = li others being zero. The matrix 
representation of the general functional (15) and the relations that must exist between the 


coefficients are studied in §5.1. 


4.2 Hybrid Forcing Potentials 

Variational principles of linear elasticity are constructed by combining the volume in- 
tegral (15) with the forcing potential P. Two forms of the forcing potential, called P d 
an< j p l jn the sequel, are of interest in the hybrid treatment of interface discontinuities. 
The d-generalized (displacement-generalized) forcing potential introduces an independent 
boundary displacement field d over S,: 

P ,l { u,<r,d) = (b,u)v + — d]s* + + [< 7 , 1 , u — d]s t . (16) 

The 1 -generalized (traction generalized) forcing potential introduces an independently var- 
ied traction displacement field t over £,■: 


P‘(u,a,t) = (b,u)v + [t,u — d)s rf + (t,u] 5 t + [t,u]s.. (17) 
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The “conventional” form P° of the forcing potential is obtained if the interface integral 
vanishes and one sets (t = <r M ] s . If so P* and P d coalesce into P c , which retains only two 
independent fields: 

P c (u,<r) = (b,u)y + ( 18 ) 


4.8 Modified Forcing Potentials 

Through various manipulations and assumptions detailed in (10] the forcing potential P 
may be transformed to 

P rf (u,ir,d) = (b,u)y + (t,d] St + (<r n ,u - d] 5 . (19) 

where the all-important surface dislocation integral is taken over S rather than S',-. One of 
the assumptions is that displacement boundary conditions (6) are exactly satisfied. This 
expression of P d is used in the sequel. A similar technique can be used to modify P*, but 
that expression will not be required in what follows. 

4.4 Complete Functionals 

Complete elasticity functionals are obtained by combining the generalized strain energy 
with one of the forcing potentials. For example, the d and t generalized versions of the 
Hu-Washizu functional are 

K = u w - P d , n t w = u w - p*. (20) 

where Uw is obtained by setting J 22 = J 13 = li J 12 = — 1> others zero, in (15). 


5. MATRIX REPRESENTATION OF ELASTICITY FUNCTIONALS 


The generalized strain energy (15) can be presented in matrix form as* 


r 

in 

J 12 

3 13 

feM 

U = A (b o c o u ) 


322 

3 23 

< e ) 

Jv 

symm 


3 '33 _ 

[e u J 


( 21 ) 


The symmetric matrix 


J 11 

312 

3 13 

symm 

322 

3 23 
3 33 


( 22 ) 


characterizes the volume portion of the variational principle. Using the relations o e = Ee, 
= EDu, e* = E -l e, and e M = Du, the above integral may be rewritten in terms of the 
independent fields as 


r 

J11E" 1 

J 12 I 

J13D 

( ~ a \ 

11 

iOk- 

3 12I 

322 E 

J23ED 

{ e ) dV. 

Jv 

.ii 3 D T 

3 23 D T E 

J33D T ED _ 

UJ 


(23) 


* To justify the symmetry of J note, for example, that j\s{d t e*)v = -3 + 

,o u )v, and so on. 



5.1 First Variation of Generalized Strain Energy 
The first variation of the volume term (15) may be presented as 

6U = (Ae,6rr) v + (Aa,6e)v - (div<r',$u)v + [<r' n ,6u]s. (24) 


where 


Ae = ji i e a + j'i 2 e + jisc”, 
A<r = ji 20 + jtt<r e + j23« rU » 
o' = ji$d + j23° e + J33 ff “- 


(25) 


The last two terras combine with contributions from the variation of P. For example, if 
p = p c the complete variation of Il c = U — P c is 


tfll c = (Ae,£ir)v + (A<r,6e)v - (div <r' + b, $u)v + -(« d,6ir n )s*. (26) 

Using P d or P* does not change the volume terms. The Euler equations corresponding to 
P d and P l are studied in (10,11) for a more restrictive form of functionals U. 

Since the Euler equations associated with the first two terms are A<r = 0 and 
Ae = 0, these quantities may be regarded as deviations from stress-balance and strain- 
compatibility, respectively. For consistency of the Euler equations with the field equations 
of §2 we must have Ae = 0, A«r = 0 and o' = o if the assumed stress and strain fields 
reduce to the exact ones. Consequently 


Jll + J 12 + Ji3 = 

Jl2 + J22 + j23 = 0, (27) 

Jl3 + J23 + J33 = 1- 

Because of these constraints, the maximum number of independent parameters that define 
the entries of J is three. 


5.2 Specific Functionals 

Expressions of J for some classical and parametrized variational principles of elasticity 
are tabulated below. The subscript of J is used the identify the functionals, which are 
listed roughly in order of ascending complexity. The fields included in parentheses after 
the functional name are those subject to independent variations. 


Potential energy (u): 


Jr 


0 0 o' 

0 0 0 . 

0 0 1 


Stress-displacement Reissner, also called Hellinger-Rcissner, (ff,u): 


J n = 


-1 

0 

1 


0 

0 

0 


1 

0 . 
0 


(28) 


(29) 
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Unnamed stress-displacement functional listed in Oden and Reddy [13 j (£,u): 




10-1 
0 0 0 . 
-10 2 


Strain-displacement Reissner-type [13) (e, u): 


( 30 ) 


Hu-Washizu (<r,e,u): 


J 5 


0 0 0 

0-11. 
0 10 


Jw 


0-11 
-1 10. 
1 0 0 


(31) 


(32) 


One-parameter stress-displacement family (<r, u) that includes Up , 
cases [9,10,11]: 


J,= 


-If 0 7 

0 0 0 
7 0 1-7 


Ur and Uu as special 


(33) 


One-parameter strain-displacement family (e, u) 

[»]: 




0 0 

0 -0 

0 0 


that includes U p and Us as special cases 


0 

0 

1-0 


(34) 


Two-parameter strain-displacement family (<r,e,u) that includes Up and t / 7 as special 
cases [9]: 

J/9-r = (1 - 0)3 , + (1 - 7 )J p -{1-0- 7)J r 


-7(1 ~0) 0 7(1 - 0) 

0 -0{1 - 7 ) 0{1 ~ 7 ) 

7(1-/?) >9(1 - 7) 1-0-1 + 20m 


(35) 


Three-parameter (a, 0, 7 ) family (a,e,u) that includes Uw and Up 7 as special cases [9]: 


= a3w + (l — <*)J/J<r 

-7(1 - 0){1 - a) —a 

= -a a- 0(1- 7)(1 - a) 

a + 7(1 -/?)(! - a) ^(1 — 7) (1 — n ) 


ot + 7 (l -/3)(1 - a) 
0(1 — 7 ) (l - a) 
(l-0-l + 20i)(l-a) 


(36) 


The last form, which contains three independent parameters, supplies all matrices J that 
satisfy the constraints ( 21 ). It yields stress-displacement functionals for a = 0 = 0 , strain 
displacement functionals for 0 = 7 = 0, and 3-field functionals otherwise. A graphic 
representation of J a p n in (a,/?, 7 ) space is given in Figure 3. 
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5.8 Energy Balancing 

A prime motivation for introducing the j coefficients as free parameters is optimization of 
finite element performance. The determination of “best” parameters for specific elements 
relies on the concept of energy balance. Let U{e) = J(E£,£)v denote the strain energy 
associated with the strain field e. If E is positive definite, U(t) is nonnegative. We may 
decompose the generalized strain energy into the following sum of strain energies: 


U = j 33 Z/(e“) + ciU{e a - e) + c 2 U{e - e“) + c 3 U{e u - e"), (37) 

where Uplc u ) = Up is the usual strain energy, c t = ^(in + J 22 - J 33 + l)» c 2 = + 

j 22 + jaa - 1), and c 3 = |(jn - J 22 + J 33 - 1). Equation (37) is equivalent to decomposing 

J into the sum of four rank-one matrices: 


J = J 33 


0 

0 

0 


0 

0 

0 


0 

0 

1 


+ Cl 


1 

-1 

0 


-1 

1 

0 


0 

0 

0 


+ c 2 


0 

0 

0 


0 

1 

-1 


0 

-1 

1 


+ C 3 


1 

0 

-1 


0 

0 

0 


-1 

0 

1 


. (38) 


Decompositions of this nature can be used to derive energy balanced finite elements by 
considering element “patches” under simple load systems. This technique is discussed for 
the one-parameter functionals generated by (34) in [5,7,8]. 
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6. FINITE ELEMENT DISCRETIZATION 

In this section assumptions invoked in the finite element discretization of the functional IT* 1 
for arbitrary J are stated. Following usual practice in finite element work, the components 
of stresses and strains are arranged as one-dimensional arrays whereas the elastic moduli in 
E are arranged as a square symmetric matrix. In the sequel we shall consider an individual 
element of volume V and surface S : Si U S t i U where Si is the portion of the boundary 
in common with other elements. 

6.1 Boundary Displacement Assumption 
The boundary displacement assumption is 

[d = Na v]s- (39) 

Here matrix Nj collects the boundary shape functions for the boundary displacement d 
whereas vector v collects the degrees of freedom of the element, also called the connectors. 
These boundary displacements must be unique on common element boundaries. This 
condition is verified if the displacement of the common boundary portion is uniquely 
specified by degrees of freedom located on that boundary. There are no derived fields 
associated with d. 

6.2 Internal Displacement Assumption 

The displacement assumption in the interior of the element is 

(u = N M q)v„ (40) 

where matrix N a collects the internal displacement shape functions and vector q collects 
generalized coordinates for the internal displacements. The assumed u need not be con- 
tinuous across interelement boundaries. 

The displacement derived fields are 

(e u = DNq = Bq) v , (<r M = EBq) v . (41) 

To link up with the FF and ANS formulations, we proceed to break up the internal 
displacement field as follows. The assumed u is decomposed into rigid body, constant 
strain, and higher order displacements: 

u = N r q r + N c q c + N,.q ; ,. (42) 

Applying the strain operator D = ^(V + V T ) to u we get the associated strain field: 

e M = DN r q r + DN c q c + DN fc q fc = B r q r + B c q c + B;,q h . (43) 

But B r = DN r vanishes because N r contains only rigid-body modes. We are also free to 
select B c = DN C to be the identity matrix I if the generalized coordinates q c are identified 
with the mean (volume-averaged) strain values e u . Consequently (44) simplifies to 

e M = e u + eJJ = e“ + B;,q fc , (44) 
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in which 


q c = 6“ = (e t ‘)v/u, (B fc )v = 0. (45) 

where v = (lW is the element volume measure. The second relation is obtained by 
integrating (44) over V and noting that q h is arbitrary. It says that the mean value of the 
higher-order displacement-derived strains is zero over the element. 

6.S Stress Assumption 

The stress field will be assumed to be constant over the element: 

{b = o)v ( 46 ) 

This assumption is sufficient to construct high-performance elements based on the free 
formulation [1-10]. Higher order stress variations are computationally effective if they are 
divergence free [10] but such a requirement makes extension to geometrically nonlinear 
problems difficult. The only derived field is 

(e ,7 =E ~ x a) v ( 47 ) 


6.4 Strain Assumptions 

The assumed strain field e is decomposed into a mean constant strain § and a higher order 

variation: > t Aa \ 

fe = e + AaV. I 48 ) 


(e = e + Aa)y. 

where e = (eV/o, A collects higher order strain modes with mean zero value over the 

e,emCnt: (A)y=0, (*>) 

and a collects the corresponding strain parameters. The only derived field is 

(<r e = Ee = Ee + EAa) v . t 50 ) 

7. UNCONSTRAINED FINITE ELEMENT EQUATIONS 

For simplicity we shall assume that all elastic moduli in E are constant over the element. 
Inserting the above assumptions into 11“ with the forcing potent.al(19) we obtama 
quadratic algebraic form, which is fairly sparse on account of the conditions (45) and (49). 
Making this form stationary yields the finite element equations 


jiivE 

juvl 

0 

-ft 

jisvl - Pi 

-p? 

j l7 vl 

jaatiE 

0 

0 

j as vl 

0 

0 

0 

jt aCfc 

0 

0 

J2sR 

— P r 

0 

0 

0 

0 

0 

■ ls vl — P« 

0 

0 

0 

jasvE 

0 

-Pa 

0 

jisR 

0 

0 


L 

0 

0 

0 

0 

0 
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where 


K, fc = (BlEBfc)v = Kj k> C,,= (A t EAV = C t , R = (B^EA) v , 

L = [NL] s , P r = [Nf n ]s, Po = [N ^ n }s, P* = [NL Is, (52) 

f r = (N^b)v, f fl = (Njb) v , fh = (N^b)v, f„ = [Njt } St . 

in which N,^ denotes the projection of shape functions N4 on the exterior normal n, and 
similarly for N r , N c and N k . Coefficient matrix entries that do not depend on the j's 
come from the last boundary term in (19). 


7.1 The P matrices 

Application of the divergence theorem to the work of the mean stress on e u yields 

(a,e u )i' = (a, e u + B^q/Jv = vd T e u + tf T (B k )vq k = vo T e u 

= [9 nt u] s = [a n ,N r q r + N c e“ + N h q h ] s = * T (P r q r + P c e“ + T h q h ). 

Hence P r = 0, P c = vl, P h = 0, and the element equations simplify to 


(53) 


ji 1 v E 1 

j'ijoI 

0 

0 

(;‘i 8 - I)vl 

0 



r a ' 


r 0 ' 


J22 

0 

0 

Jaaul 

0 

0 


e 


0 

0 

0 

J22C* 

0 

0 

j*R T 

0 


a 


0 

0 

0 

0 

0 

0 

0 

0 

4 

q r 

► = i 

fqr > 

(jis - l)vl 

J 2 SVI 

0 

0 

jaauE 

0 

0 


e u 


f-Jtt 

0 

0 

J23R 

0 

0 

J33K q h 

0 


<lh 


Igh 

L 

0 

0 

0 

0 

0 

0 _ 

1 

< v , 


5 Tv j 


( 54 ) 


The simplicity of the P matrices comes from the mean-plus-deviator expression (44) for 
e u . If this decomposition is not enforced, P r = 0 but P c = (B c )v and P^ = (B/Jv 

8. KINEMATIC CONSTRAINTS 

The “tricks” we shall consider here are kinematic constraints that play a key role in the 
development of high-performance FF and ANS elements. These are matrix relations be- 
tween kinematic quantities that are established independently of the varijitional equations. 
Two types of relations will be studied. 

8.1 Constraints Between Internal and Boundary Displacements 

Relations linking the generalized coordinates q and the nodal connectors v were introduced 
by Bergan and coworkers in conjunction with the free formulation (FF) of finite elements 
(2-3). For simplicity we shall assume that the number of freedoms in v and q is the same; 
removal of this restriction is discussed in [10]. By collocation of u at the element node 
points one easily establishes the relation 

v = G r q r + G c q c + G k q k = Gq, (55) 
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where G is a. square transformation matrix that will be assumed to be nonsingular. On 
inverting this relation we obtain 


q = G" 1 = Hv, 




qr 1 


H r 

or 

q = \ 

e« 1 

► = 

H c 


1 

Ufc J 


H/i 


(56) 


The following relations between L and the above submatrices hold as a consequence of the 
individual element test performed in §9.3: 


L r G r = 0, L t G c = vl, vH c = L t . (57) 

If the decomposition (44) is not enforced, the last two should read L T G C = vB c , a relation 
first stated in [3], and P C H C + P/iH/» = L r . 

8.2 Constraints Between Assumed Higher Order Strains and Boundary Displacements 

Constraints linking e h to v are of fundamental importance in the assumed natural strain 
(ANS) formulation. The effect of these constraints in a variational framework is analyzed 
in some detail in [11-12]. Here we shall simply postulate the following relation between 
higher order strains and nodal displacements: 

a = Qv. (58) 


where Q is generally a rectangular matrix determined by collocation and/or interpolation. 
The individual element test in §9.3 requires that Q be orthogonal to G r and G c : 

QG r = 0, QG C = 0. (59) 

The constraint (58) still leaves the independently varied mean strain e to be determined 
variationally. 


9. VISIBLE STIFFNESS EQUATIONS 

Enforcing the constraints a = Qv, q r = H r v, q c = H c v = t» _l L r v, q^ = H^v, through 
Lagrange multiplier vectors A ( „ A r , A c , and A,„ respectively, we get the augmented finite 

element equations 


ji i tiE 

jlivl 

0 

0 

(iis - i)ul 

0 

0 

0 

0 

0 

L r 1 


a ' 


' 0 1 

j iaul 


0 

0 

Jasol 

0 

0 

0 

0 

0 

0 


e 


0 

0 

0 

Jaac/k 

0 

0 

J 2 S R T 

-I 

0 

0 

0 

0 


a 


0 

0 

0 

0 

0 

0 
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Condensation of all degrees of freedom except v yields the visible * element stiffness equa- 
tions 

Kv = (Kb + K/,)v = f (61) 

where 

Kb = v - 1 LEL t , (62) 

K h = j 33 Hj[K v bHb + j '23 (H^RQ + Q r R T Hb) + j 22 Q T CbQ, (63) 

f = f„ + H T r f qr + v~ l lJi qc + Hlf qh . (64) 

Adopting the nomenclature of the free formulation (3], we shall call Kj, the basic stiffness 
matrix and Kb the higher order stiffness matrix. 

9.1 Relation to Previous HP Element Formulations 

If J = J 7 of (33), y 33 = l — '■y, J 22 = y 23 = o, and we recover the scaled free formulation 
stiffness equations studied in [5,7,9,10]: 

Kb = (l-7)HjK ff bHb. (65) 

If we take J = Jw of (32), j '22 = 1, j 33 = j 23 = 0 and we obtain 

Kb = Q T CbQ. (66) 

This is similar to the stiffness produced by the ANS hybrid variational formulation studied 
in [11-12], in which the potential P l was used instead of P d . 

But the term with coefficient j 2 3 in (63) is new. It may be viewed as coupling the FF 
and ANS formulations. It is not known at this time whether (61-64) represents the most 
general structure of the visible stiffness equations of HP elements. 

9.2 Recovery of Element Fields 

For simplicity suppose that the body forces vanish and so do f flr , f qc and f q h.. If v is known 
following a finite element solution of the assembled system, solving the equations (60) for 
the internal degrees of freedom yields 

e = t» _1 L T v, fr = Ee, a = Qv, q r = H r v, 6“ = 6, q fc = HbV, 

A a = (j 22 C!bQ + j 33 R r Hb)v, A r = 0, A c = 0, A;» = (j 23 RQ + j33K y bHb)v. 

It is seen that the mean strains e, e u and E -1 o agree, and so would the mean 

stresses. This is not the case, however, if the body forces are not zero. It is also worthwhile 
to mention that a nonzero Lagrange multiplier vector flags a deviation of the associated 
fields from the variationally consistent fields that would result on using the unconstrained 
FE equations (54) without “tricks”. 


* The qualifier visible emphasizes that these are the stiffness equations other elements “see”, 
and consequently are the only ones that matter insofar as computer implementation on a 
displacement-based finite element program. 
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9.S The Individual Element Test 

To conclude the paper, we investigate the conditions under which HP elements based on 
the foregoing general formulation pass the individual element test of Bergan and Hanssen 
[1-3]. To carry out the test, assume that the “free floating” element* under zero body 
forces is in a constant stress state a 0 , which of course is also the mean stress. Insert the 
following data in the left-hand side vector of (60): 


* = *<) = &“, e = E"V 0 , a/v = 0, q r = arbitrary, e tt = e u = E l 9 0 , q h -0, 

_ _ u i /""* MS' “ 1 * 


= 0, A r = 0, A c = 0, A/, = 0, v = G r q r + G c e u — G r q r + G C E a 0 . 


( 68 ) 


Premultiply by the coefficient matrix, and demand that all terms on the right-hand side 
vanish but for f„ = L<r 0 . Then the orthogonality conditions in (57) and (59) emerge. This 
form of the patch test is very strong, and it may well be that relaxing circumstances can 
be found for specific problems such as shells. 


10. CONCLUSIONS 

The results of the present paper may be summarized as follows. 

1. The classical variational principles of linear elasticity may be embedded in a 
parametrized matrix form. 

2. The elasticity principles with assumed displacements are members of a three- 
parameter family. 

3. Finite element assumptions for constructing high-performance elements may be con- 
veniently investigated on this family. 

4. Kinematic constraints established outside the realm of the variational principle may 
be incorporated through Lagrange multiplier adjunction. 

5. The FF and ANS methods for constructing HP finite elements may be presented 
within this variational setting. In addition, combined forms emerge naturally from 
the general parametrized principle. 

6. The satisfaction of the individual element test yields various orthogonality conditions 
that the kinematic constraints should satisfy a prton. 

The construction of high performance elements based on a weighted mix of FF and ANS 

“ingredients” will be examined in sequel papers, and specific examples given to convey the 

power and flexibility of the present methods. 
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